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The mechanism of thermal leptogenesis is investigated in the high-energy regime of braneworld 
cosmology. Within the simplest seesaw framework with hierarchical heavy Majorana neutrinos, 
we study the implications of the modified Friedmann equation on the realization of this mecha- 
nism. In contrast with the usual leptogenesis scenario of standard cosmology, where low-energy 
neutrino data favors a mildly strong washout regime, we find that leptogenesis in the braneworld 
regime is successfully realized in a weak washout regime. Furthermore, a quasi-degenerate light 
neutrino mass spectrum is found to be compatible with this scenario. For an initially vanishing 
heavy Majorana neutrino abundance, thermal leptogenesis in the brane requires the decaying heavy 
Majorana neutrino mass to be Mi > 10^" GeV and the fundamental five-dimensional gravity scale 
10^^ < AI5 < 10^® GeV, which corresponds to a transition from brane to standard cosmology at 
temperatures 10** < Tt < lO^"* GeV. 
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I. INTRODUCTION 

The most recent Wilkinson Microwave Anisotropy 
Probe (WMAP) results and big bang nucleosynthesis 
analysis of the primordial deuterium abundance imply 1] 

= ~ = (6.1 ± 0.3) X 10-1" , (1) 

for the baryon-to-photon ratio of number densities. 
Among the viable mechanisms to explain such a matter- 
antimatter asymmetry in the universe, leptogenesis 
has become one of the most compelling scenarios due 
to its simplicity and close connection with low-energy 
neutrino physics. In the simplest framework, consist- 
ing on the addition of right-handed neutrinos to the 
standard model, leptogenesis can be easily realized by 
means of the out-of-equilibrium decays of the heavy Ma- 
jorana neutrinos Ni at temperatures below their mass 
scale Mi. The lepton asymmetry generated in the pres- 
ence of CP- violating processes is then partially converted 
into a baryon asymmetry by the sphalerons 

Right-handed neutrinos can also provide a natural 
explanation for the smallness of the neutrino masses 
through the well-known seesaw mechanism ■ In this re- 
spect, the heavy Majorana neutrino mass spectrum turns 
out to be crucial for a successful implementation of both 
leptogenesis and the seesaw mechanism. In particular, 
the standard thermal leptogenesis scenario with hierar- 
chical heavy Majorana neutrino masses (A/i ^ M2 < 
M3) constrains the lightest heavy Majorana mass to be 
Ml > 4 X 10^ GeV ^, ^ JJ, if A^i is in thermal equihb- 
rium before it decays, or Mi > 2 x 10^ GeV for a zero 



•Electronic address; | bento@sirius . ist . utl . ptj 
^Electronic address: gonzalez@cftp.ist.utl.pt] 
■f Electronic address: ncsantos@cftp.ist.utl.ptl 



initial iVi population. Moreover, an upper bound on the 
light neutrino masses, mi < 0.12 eV ,9i, is implied by 
successful leptogenesis. 

In determining the departure from thermal equilib- 
rium, the interplay between the expansion rate of the 
universe H and the particle reaction rates involved at 
the relevant epoch is crucial. Of course, such an inter- 
play depends on the specific properties of early cosmol- 
ogy. Braneworld cosmology (BC) has opened up the pos- 
sibility for a new phenomenology of the early universe. 
In particular, the Randall-Sundrum type II braneworld 
model ^3 has recently received much attention. In 
this model, the expansion dynamics at early epochs is 
changed by the presence in the Friedmann equation of a 
term quadratic in the energy density ll l| , 

where Mp ~ 1.22 x 10^^ GeV is the four-dimensional 
Planck mass, M5 is the five-dimensional Planck mass and 
we have set the four-dimensional cosmological constant to 
zero and assumed that inflation rapidly makes any dark 
radiation term negligible. Notice that Eq. ||2J) reduces to 
the usual Friedmann equation at sufficiently low energies, 
p ^ A, while at very high energies we have H oc p. This 
behavior has important consequences on early universe 
phenomena such as inflation |l2l | and the generation of 
the baryon asymmetry [l^ . 

In this paper, we study how the standard mechanism of 
thermal leptogenesis is affected by modifications in the 
expansion rate of the universe due to early braneworld 
cosmology. We perform a thorough analysis of this sce- 
nario by solving the relevant Boltzmann equations fol- 
lowing the approach of Refs. for the standard cos- 
mology (SC) case. 

The structure of the paper is as follows. In Section II 
we briefly present the basic formulae for the mechanism 
of thermal leptogenesis in braneworld cosmology. In sec- 
tion III we present a simple fit for the efficiency factor 
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as a function of the decay parameter, obtained by solv- 
ing numerically the Boltzmann equations. In Section IV 
we find lower bounds on the decaying heavy Majorana 
neutrino mass, the fundamental five-dimensional gravity 
scale and the transition temperature from brane to stan- 
dard cosmology. In Section V we derive the upper bounds 
on the light neutrino masses implied by leptogenesis. Fi- 
nally, in section VI we present our conclusions. 



II. THERMAL LEPTOGENESIS IN THE 
BRANE 

In its simplest standard model version, leptogenesis is 
dominated by the CP- violating interactions of the light- 
est of the heavy Majorana neutrinos. Assuming that 
right-handed neutrinos are hierarchical. Mi ^ Af2,3 , and 
that the decays of the heavier neutrinos do not influence 
the final value of the B — L asymmetry, the baryon asym- 
metry generated by Ni can be conveniently parameter- 
ized as 



VB = j A^i 



(3) 



where ^ = 28/79 is the fraction of the B — L asymmetry 
converted into a baryon asymmetry by sphaleron pro- 
cesses Hi and / = 2387/86 is the dilution factor cal- 
culated assuming standard photon production from the 
onset of leptogenesis till recombination. The amount 
of B — L asymmetry per comoving volume, Ng_j^ = 
Nb-l{T <C Ml), is obtained from the solution of the 
relevant Boltzmann equations, which in the above mini- 
mal framework can be written in the compact form 0, 



Ni 



dz 

dNB-L 

dz 



{D + S){Nr,,-N%\) , (4) 
eiD (Nn, - N'J ) - WNb-l , (5) 



where z — Mi /T , N'^^ is the equilibrium number density, 
normalized so that 7V^' (z ^ 1) = 3/4, and ei is the CP- 
asymmetry parameter in iVi decays. Different physical 
processes contribute to these equations. The terms D, S 
and W, defined through the corresponding rates Pi as 



D 



Pd 
Hz ' 



Hz 



account for decays and inverse decays (P), AP = 1 
scatterings {S) and washout processes (W), respectively. 
While decays yield the source for the generation of the 
B — L asymmetry, all other processes, including AP = 2 
processes mediated by heavy neutrinos, contribute to the 
total washout term. 

During radiation domination p = TT^g^,T^ /30 and the 
expansion rate (0) is given by 



H{z) 



H / l + (zt/z)4 



z 



l + zf 



(7) 



where 



H = H{z = 1) 



/ 47r3g, /— 7 



(8) 



and 5* = 106.75 is the effective number of relativistic 
degrees of freedom. The parameter zt is defined as 



Zt 



Ml 



TT^g^M^If 
45M| 



(9) 



where Pt is the temperature at which the transition from 
brane to standard cosmology takes place (p(Pt) — 2A), 



Tt ~ 9.8 X 10^ GeV 



M. 



1012 GeV 



3/2 



(10) 



In establishing a connection between leptogenesis and 
neutrino physics, the dependence of the different con- 
tributions given in Eq. © on the neutrino parameters 
turns out to be crucial. In particular, one can show that 
terms D and S are proportional to the so-called effective 
neutrino mass parameter 



mi 



{m'jjmD)ii 
All 



(11) 



where is the Dirac neutrino mass matrix. On the 
other hand, the washout term W contains in general two 
contributions, one which depends only on rhi and another 
which is proportional to the absolute neutrino mass scale 
= -|- 7712 + "^3 ■ 
Departure from thermal equilibrium is controlled by 
the decay parameter K , 



K 



Pp _ 7771 

H 



(12) 



where Pp = Pp(z = oo) is the iVi decay width at zero 
temperature, 



P, 



D 



1 7771 



(13) 



87r ^2 ' 

u w 174 GeV and 777* is the equilibrium neutrino mass 

777* = ^T^H^ ^ 1.08 X 10-^ eY^Jl + zf. (14) 

In the high-energy regime of brane cosmology, i.e. for 
Zt ^ 1, we find 



1.1 X 10"^ eV 



Ml 



^ ^10^2 GeV Y 



IQSGeVy V J 
The solution of Eq. Q for Nb-l is given by 



(15) 



Nb-l{z) = N^B-Le--^'^'''"^^''^ - \ ei n{z) , (16) 
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where Ng_j^ accounts for any possible initial B — L asym- 
metry (e.g. due to the decays of the heavier neutrinos 
■^2,3) and k{z) is the efficiency factor, which measures 
B — L production from iVi decays, 

Assuming no pre-existing asymmetry, i.e. Ng^i^ = 0, 
from Eqs. lj3Jl and (|16|l the final baryon asymmetry can 
be conveniently written as 

rjB ^ dtiKf , (18) 

where d ~ 0.96 x 10^^ and k/ = k(oo) is the final value 
of the efficiency factor, normalized in such a way that 
K/ ^ 1 in the limit of thermal initial iVi abundance and 
no washout. The computation of is a rather difficult 
task, since it involves solving the Boltzmann equations. 
Nevertheless, simple analytical formulae can be found in 
the limiting cases of strong {K 1) and weak {K -^1) 
washout regimes 0, and numerical fits valid for all val- 
ues of ffu can also be derived in the standard cosmology 
case 13, 1^ . Clearly, one expects these results to be mod- 
ified in the braneworld scenario. In the next section we 
will derive a simple fit for the efficiency factor, assuming 
that the baryon asymmetry is created during the high- 
energy regime of brane cosmology. 

III. THE EFFICIENCY FACTOR 

Neutrino oscillation experiments presently constrain 
two mass squared differences for the light neutrinos [l5j |. 



nisoi = ^J l^ml^i ~ 8.9 x 10"=^ eV , 

mat™ = ^/a^^~ 5.0 X 10-2 eV, (19) 

from solar and atmospheric measurements, respectively. 
The most stringent bounds on the absolute neutrino mass 
scale come from Oz//3/3-decay experiments |l(|, which 
yield an upper bound on the light Majorana neutrino 
masses of about 1 eV, and from the 2df galaxy red- 
shift survey and WMAP results Q, which indicate that 
< 0.7 eV. 

The effective neutrino mass parameter mi is an essen- 
tial quantity for thermal leptogenesis. Remarkably, it 
is possible to find a natural range for this quantity in 
terms of the light neutrino masses. Indeed, the lower 
bound m-i > mi is always verified |l7j |. Moreover, in 
the absence of strong cancelations or specific fine-tunings, 
one has the upper bound mi < Q. Thus, one finds 
mi < mi < m^. 

Of particular interest is the favored neutrino mass 
range nisoi ^ wi ^ ruatm- In this situation, leptogen- 
esis in the standard cosmology case would occur in the 
mildly strong washout regime 8 < K < A6 (see Eqs. (|12f) 
and II4I1 I. On the other hand, in the high energy regime 
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FIG. 1: (Color online) The evolution of Nni (red curves), A''^''^ 
(green curve) and |7?B/ei| (black curves) as functions of the 
temperature for Mi = 10^^ GeV and fhi = niatm- The solid 
(dashed) hues correspond to the brane (standard) cosmology 
result. 

of brane cosmology {zt ^ 1) one has 8/ zf < K < 46/ , 
which for zt > 7 already implies that leptogenesis in 
the brane lies in the weak washout regime. As an il- 
lustrative example, in Fig. ^ we compare the evolution 
of the abundances with temperature in the SC strong 
and BC weak washout regimes. Assuming zero initial 
A^i population at T ^ Mi, we have numerically solved 
the system of Boltzmann equations jl)) and (O follow- 
ing the approach presented in Ref. [g, which properly 
includes the leading finite-temperature corrections and 
takes into account all the relevant processes. The evo- 
lution of the A^i abundance (red curves) and the ratio 
l^ys/cil (black curves) are shown as functions of z for 
Ml = 10" GeV and mi = m^t™ ~ 0.05 eV. The dashed 
lines correspond to the standard cosmology result, while 
the solid lines are the brane cosmology result, taking 
M5 = 4.5 X 10^2 GeV. With this choice of parameters, 
Ksc = 46 and Kbc = 4 x 10"^. 

It is possible to find a simple numerical fit for the effi- 
ciency factor K/ in terms of the decay parameter K. In 
both SC and BC cases, the behavior of k/ is well de- 
scribed by the power law fit 

valid for a heavy neutrino mass Mi < 10^** GeV. For SC 
we find 

a = 3.5, 6 = 0.6, c= 1.2, (21) 

(see also Ref. 83), while in the high-energy regime of BC 
we obtain 

a = 0.8, 6 = 0.025, c= 1.65. (22) 
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FIG. 2: (Color online) The efficiency factor as a function 
of the decay parameter obtained by solving numerically 
the Boltzmann equations, for different values of the effective 
neutrino mass mi. The blue and red solid curves correspond 
to the numerical fits of «:/ given in Eqs. (Iinj, lHU and ll^ . 
for the cases of standard and brane cosmology, respectively. 



These fits are plotted in Fig. [3 and compared with the 
efficiency factor obtained by solving numerically the full 
set of Boltzmann equations for different values of rh\. As 
expected, the curves in the BC case are displaced to the 
weak washout region where thermal leptogenesis is more 
efficiently realized. 

Using the numerical fit (|20|l . one can see that maximal 
efficiency is achieved at 



-^peak — 



c 



c\ l/(c+l) 



which in BC corresponds to 

Kpeak ^ 0.08 , nf^"" ~ 0.06 , 
while in SC one finds 



peak 



1.23, 



peak 

V 



0.19. 



(23) 



(24) 



(25) 



We notice that this value of in the standard cos- 
mology case implies an effective neutrino mass fhi ~ 
1.3 X 10~^ eV. On the other hand, the maximal k^^^ 
in the high-energy regime of BC can be achieved for any 
value of fhi > 10""^ eV, including the favored neutrino 
mass range rusoi < mi < matm- 

The above values for the maximal efficiency factor to- 
gether with the upper bound on the CP asymmetry ei 
can be translated into bounds on the heavy Majorana 
neutrino mass Mi , the 5D Planck scale M5 , the transi- 
tion temperature Tj and the light neutrino masses rrii . 



IV. LOWER BOUNDS ON Mi , M5 AND Tt 

In order to discuss the bounds implied by leptogenesis, 
it is convenient to write the CP asymmetry ei as the 
product 



ei 



max 
^1 



sin St 



(26) 



where e™'^'^ is the maximal CP asymmetry and 6l is an 
effective leptogenesis phase. In general^e™'''' depends on 
Ml, fhi and the light neutrino masses [g- In the case of 
hierarchical and quasi-degenerate light neutrinos. 



3 Ml 
167r 



1 




(27) 



where for a normal hierarchy 



TO, 



— 2 o 2 

TO =6 171, 



'•sol J 



(28) 



2 m 



sol ' 



The CP-asymmetry reaches its maximum for mi = 
and ma ~ niatm j i-e- for fully hierarchical neutrinos. In 
this case, 



3 MirUa 
167r 

10"^ 



Ml 



IQi" GcV 



/ mat 
V0.05 



eV 



(29) 



A lower bound on Mi can be obtained from the con- 
dition 



max 

Vb 



= del 



max ^pcak 



^ Vb J 



(30) 



where 77^'''' is given in Eq. and rj'^^^ is the maximal 
baryon asymmetry with the efficiency factor k^"^*^ given 
in Eqs. H24|l and H25|l for the BC and SC regimes, respec- 
tively. 

The condition l|30(l then leads to the following lower 
bound on Mi 



Ml > 6.3 X 10® GeV 



obs 

Vb 



6 X 10 



-10 



0.05 eV 



(31) 



Taking into account the minimum value for r]B allowed 
by observations 77^'''' = 5.8 x 10"^" (cf. Eq. (P)) and 
the atmospheric neutrino mass difference result given in 
Eq. H19() . one obtains 



Ml > 6.1 X 10^ GeY {kY'-^)-\ 



(32) 



In the SC leptogenesis regime, k^^^ — 0.19 and one gets 
Ml > 3.2 X 10^ GeV. On the other hand, in the high- 
energy regime of brane cosmology, k^'^^^ ~ 0.06 and one 
obtains the more restrictive bound 



Ml > 1.0 X 10'"' GeV. 



(33) 
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This bound implies, in turn, a lower bound on the 
gravity scale M5. Indeed, imposing K — i^Tpoakj with K 
given in Eq. (fT^ . and using Eq. and the bound (p?^ . 
we get 



M5 > 9.7 X lO^-^ GeV 



mi 



10-3 eV 



-1/3 



(34) 



Finally, from this result and using Eq. (|10|l . we can derive 
the following lower bound on the transition temperature 



> 3.0 X 10^ GeV 



mi 



-1/2 



10-3 eV, 

In the favored neutrino mass range msoi ^ 



(35) 



mi < TTlatm, 



the bounds (|34|l and (|35() imply 

A/5 > (2.6 - 4.7) X 10^2 , (36) 
Tt > (4.2 - 9.9) X 10^ GeV . (37) 

For mi < 1 eV, these bounds yield A/5 > 1.0 x lO^^ Qgy 
and Tt > 1.0 X 10^ GeV. 
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FIG. 3; (Color online) Contour lines of constant r?a'"' = 
rj°j^'' = 5.8 X 10-^° in the (mi , M5)-plane for different val- 
ues of mi . We assume that thermal leptogenesis occurs in 
the high-energy regime of brane cosmology. 



V. UPPER BOUND ON THE LIGHT 
NEUTRINO MASSES 

The washout term W receives contributions from in- 
verse decays, AL = 1 processes and AL — 2 pro- 
cesses. At low temperatures {z 3> 1) the contribution 
AW, which depends on the absolute neutrino mass scale 
m, dominates over the contribution proportional to m-i. 
This contribution is approximately given by 



AW{z > 1) 

where 

9\/5A/pA/im2 



cw 



(38) 



Cw 



8 7r9/2^w4 



0.19 



A/i 



1010 GeV / V 1 eV 



X 2 

m 



(39) 



The presence of AW modifies the total efficiency fac- 
tor. 



_ - f °° dzAW 

Kf = Kf e 'B 



(40) 



where zb is the value at which the lepton asymmetry is 
no longer produced. The washout term in Eq. H38I) leads 
then to 



Kf = Kfe 



-Cw f(zB,Zt) 



with 



f{zB,Zt) 



dz ■ 



1 

z^^l + {zt/z) 

4 



B 



zt 



1 1 

'4' 2 



(41) 

(42) 
(43) 



where /?(z; a, 6) is the incomplete beta function. For zt <C 
Zb one recovers the standard cosmology result (2j 



!{zb)^—, 
Zb 

while for zt 3> (BC regime) one finds 
r(l/4)r(5/4) 1.85 



!{zt) 



/TT Zt 



Zt 



(44) 



(45) 



leading to 



Kf ~ Kf exp 



-3.4 X 10" 



Mr, 



1012 GeV 



3/2 / 



leV 



(46) 



We notice that the dependence of the washout term on 
A/i is canceled out in the BC regime. In the latter case, 
the AW contribution to the total washout depends only 
on the fundamental scale A/5 and the absolute neutrino 
mass scale m. Clearly, the exponential suppression in the 
efficiency factor becomes relevant only for large values of 
A/5 or m. In Fig.|3|we present the contour lines of con- 
stant 7?^^'= = ■q°j^'' = 5.8 X 10-1° in the (mi , A/5)-plane 
for different values of m-i . The baryon asymmetry if^^^ 
was computed using the maximal CP-asymmetry (cf. 
Eq. 123) and taking into account the efficiency factor k/ 
as given in Eq. (|46|l with k f at its peak value H24|) . From 
the figure one finds the upper bound A/5 < lOi^ GeV, 
which translates into an upper bound on the transition 
temperature, Tt < lOi^ GeV. On the other hand, the 
light neutrino mass mi turns out to be constrained es- 
sentially by the value of mi. For mi < 1 eV, we find 
mi < 0.8 eV. Thus, one concludes that the thermal lep- 
togenesis scenario in braneworld cosmology is consistent 
with quasi-degenerate light neutrinos. 
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VI. CONCLUSIONS 

The explanation of the matter-antimatter asymmetry 
observed in the universe remains an open question. Lep- 
togenesis is a simple and elegant mechanism for explain- 
ing the cosmological baryon asymmetry. Since this asym- 
metry must has been created at early times, a successful 
realization of this mechanism crucially depends on the 
properties of early cosmology. 

In this paper, we have studied the implications of 
Randall- Sundrum type II braneworld cosmology on the 
thermal leptogenesis scenario. In contrast with the usual 
leptogenesis scenario of standard cosmology, where low- 
energy neutrino data favors a strong washout regime 
and quasi-degenerate light neutrinos are not compatible 
with leptogenesis, we have found that in the high-energy 
regime of braneworld cosmology leptogenesis can be suc- 
cessfully realized in a weak washout regime, and a quasi- 
degenerate neutrino mass spectrum is allowed. 

For an initially vanishing heavy Majorana neutrino 
abundance, we have obtained bounds on the decaying 
heavy Majorana neutrino mass, the fundamental five- 
dimensional gravity scale and the transition temperature 
from brane to standard cosmology. As far as the up- 
per bound on the light neutrino masses is concerned, we 
have seen that thermal leptogenesis in the brane imposes 
a limit less restrictive than rrii < 0.23 eV, which is the 
one presently implied by WMAP results. Finally, it is 
worth emphasizing that all the bounds have been de- 
rived under the simple assumption of hierarchical heavy 
Majorana neutrinos. For a partially degenerate spectrum 
these bounds can be relaxed. Indeed, if at least two of 
the heavy Majorana neutrinos are quasi-degenerate in 
mass, i.e. Mi « M2, then the leptonic CP asymmetry 
relevant for leptogenesis exhibits the resonant behavior 
ei - Mi/{M2 - Ml) [3, and values of ef'''' - can 
be reached. In this case, the lower bound on Mi coming 
from leptogenesis can be relaxed H, . 

Let us also briefly comment on the consistency of infla- 
tion and reheating with the present framework of thermal 
leptogenesis. Assuming that inflation is driven by a sim- 
ple quadratic chaotic inflation potential, V = ^m^ip'^, the 
scale of inflation in the p^-dominated period is roughly 
given hy V ^ (0.1 Ms)^ thus ensuring that sufficient 
inflation can take place while fulfilling the requirement 
V < . Another important constraint comes from the 
reheating process, often associated with the decay of the 
inflaton into elementary particles, and which eventually 
results in the creation of a thermal bath. If the inflaton 
decays through Yukawa interactions with a coupling to 
matter fields of order ~ 0{1), then the decay rate is fast, 
~ TO, and reheating is practically instantaneous. As- 
suming that inflation ends in the p^-regime, the reheating 
temperature is estimated as Trh ~ 0(10"^) x M5 , so that 
the energy density of radiation prh = Tr'^g^T^i^/SO < M| . 



Since consistency with braneworld thermal leptogenesis 
requires Tt < Mi < Trh , from Eq. I|10(l one obtains the 
upper bound A/5 < 10^^ GeV, which coincides with the 
leptogenesis bound implied by the observed baryon asym- 
metry (cf. Fig. 3). Notice however that, for values of M5 
close to the upper limit, i.e. M5 ^ 10^^ GeV, one would 
have Tt ~ Mi ~ Trh — lO^'' GeV and the mechanism 
of leptogenesis would occur during the transition from 
brane to standard cosmology. 

During the reheating process the distribution of par- 
ticles is far from thermal equilibrium. Thus, number- 
conserving as well as number-violating interactions are 
necessary to achieve kinetic and chemical equilibrium 
among the different species. In an expanding universe, 
full thermal equilibrium requires a sufficiently low tem- 
perature. Indeed, 2^2 GUT interactions mediated 
by massless gauge bosons would have occurred in the 
very early universe at a rate F ~ Q^qut with acuT ^ 
0(10"^). Since during the radiation-dominated era of 
standard cosmology H ~ T'^/Mp , such interactions are 
expected to be in equilibrium for T < Q^gut ^^p ^ 
10^^ GeV J9j. On the other hand, in the high-energy 
regime of braneworld cosmology the expansion rate scales 
as ~ T^ /M^ and, consequently, thermal equilibrium 
is expected at T < a^^r^M^ < 5 x 10^^ GeV for 
M5 < 10^^ GeV. Clearly, the above analysis could be sub- 
ject to modifications, if there exist new interactions which 
could thcrmalizc the universe at the earliest epochs. For 
instance, it is conceivable that quantum gravity effects 
near the fundamental Planck scale could have resulted in 
a state of maximal entropy. Nevertheless, in the absence 
of a complete theory of quantum gravity, this remains an 
open issue. 

Note added: While this work was in preparation, a 
related preprint has appeared [2(1 |. However, the conclu- 
sions drawn by the authors about mass scales and param- 
eters somewhat differ from ours because in Ref. [23 the 
corresponding bounds were derived by using the right- 
handed neutrino thermalization condition instead of solv- 
ing the full set of Boltzmann equations. 
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